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Applying standard mathematical methods, it is explicitly shown how the Riccati equation 
for the Hubble parameter H(rj) of barotropic open FRW cosmologies is connected with 
a Kortcwcg-de Vrics equation for adiabatic index solitons. It is also shown how one can 
embed a discrete sequence of adiabatic indices of the type n 2 (|7 — l) 2 (7 7^ 2/3) in the 
sech FRW adiabatic index soliton. 



1 - Riccati equation for barotropic FRW cosmologies 

Recently, FaraoniJ3 showed that the equations describing barotropic FRW cosmolo- 
gies can be combined in a simple Riccati equation leading in a much easier way to 
the textbook solutions for the FRW scale factors. Faraoni obtained the following 
cosmological Riccati equation 

— — = —cH — kc , (la) 

for the log derivative of the FRW scale factor, the famous Hubble parameter H(rf) = 
da l^ , The independent variable is the conformal time 77, c is simply related to the 
adiabatic index of the cosmological fluid, c = §7 — 1 which was assumed constant 
by Faraoni, and K = 0, ±1 is the curvature index of the flat, closed, open FRW 
universe, respectively. For the mathematical considerations of this work we are 
more interested in the open cosmologies n = — 1. The flat case is special and has 
been discussed in a previous work^ whereas the case K = +1 will be commented in 
subsection 3.2 herein. 

We notice that from the mathematical point of view, Faraoni's FRW Riccati 
equation being of constant coefficients is directly integrable. The solutions are 

H + = — tan(c?7 — <f> + ) , H~ = tanhfc?? — <£_), (16) 

for the closed and open FRW universes, respectively. Of course, depending on 
the initial conditions, the solutions can be written in terms of the corresponding 
cofunctions as well. The minus sign in front of the phases has been chosen for later 
convenience. 
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2 - Cosmological Korteweg de Vries equation 

We now pay special attention to the k = — 1 case for which a connection with the 
KdV equation is possible. Employing the well-known change of function H~ = 
one can pass from the nonlinear Riccati equation to the linear second-order 
differential equation 

w" - c 2 w = , (2) 
where the prime means derivation with respect to rj. This equation can be factorizcd 

A 1 A 2 w = (3) 

by means of the operators 

^-{i +w )- (4 » 

where W = cH~ . The general linear solution (the general zero mode) is of the form 

w(ij, T) = a(T)e c " + /?(!>-"' , (5) 

where we have assumed that the superposition constants are functions of a new 
coordinate T that is independent of r\. It is easy to show the relationship between 
the two superposition constants and the phase <fi- 

0_(T) = iln/3/a. (6) 

Inverting as in SUSY QM the order of application of the factorization operators, 
one gets an equation of the type 

A 2 A lU = or u" + c\{ri, T)u = Q , (7) 

where Ci(rj,T) is the sech soliton with the asymptotic wings going to c 2 

cl(r],T) = c 2 - 2c 2 sech 2 ( C 77 - </»_(T)) . (8) 

It is now easy to apply basic results from the KdV mathematics allowing to interpret 
c\ as a single sech soliton of the following KdV equationJ§ 

d C „„dC d 3 C n 

if the condition -fc = 4c 3 is imposed. Moreover, by choosing initial conditions at 
T = 0, one can see that T is the evolution variable whereas the conformal time 
plays a role similar to a spatial variable Ji 
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3 - SUSY recursive procedure 

3.1 - The scheme presented in the previous section can be iterated in a simple and 
well-known manner.™ We put T — const, i.e., we consider an instantaneous soliton 
profile and take 0_ = for simplicity. Thus, we start with a 'quantum' mechanical 
system in a constant 'potential' c 2 and relate it to a Schr6dinger-type equation in 
the conformal time domain which has a fundamental 'frequency' at — c 2 . We solve 
the "fcrmionic" Riccati equation i.e., 

W? - W[ - c 2 = , (10) 

to find Witten's superpotential W\[rj) — — ctanh[c?7] and next go to the "bosonic" 
Riccati equation 

W 2 + W[ + c 2 ^) - c 2 = , (11) 

in order to get c 2 (rj) = — 2c 2 sech 2 [cr;]. Moreover, one can write the Schrodinger 
equation corresponding to the "bosonic" Riccati equation as follows 

-y" + c\{ri)y = -c 2 y , (12) 

with the solution y oc csech(c7y). The physical picture is that of a sech soliton 
profile containing a single mode self-trapped at — c 2 within the frequency pulse. 
One can employ the scheme recursively to get several localized modes embedded in 
the soliton profile. Indeed, suppose we would like to introduce N adiabatic indices 
of the type c 2 = — n 2 c 2 , n = 1, ...N in the sech pulse. Then, one has to solve the 
sequence of equations 

W? l -Wi = c 2 n _ 1 +n 2 c 2 , (13) 
W 2 + W' n = c 2 n + n 2 c 2 , (14) 

inductively for n = l...i\r.@@ The soliton containing TV adiabatic indices n 2 c 2 , n = 
1...N is of the form c 2 N {rf) = —N(N+ l)c 2 sech 2 (c?y). The corresponding modes can 
be written in a compact form as follows 

y„(r);N) w A^(r ] ;N)A'<(ri;N-l)Al(T r ,N-2)...AHr);N-n + 2)scch N - n+1 cr) , (15) 

i.e., by applying the first-order operators A^(ri;a n ) = — -j- — a„ctanh(c?y), where 
a n = N — n, onto the "ground state" sech mode. 

3.2 - For the case k = 1 the "fermionic" Riccati equation 

W 2 - W' x + c 2 = (16) 
leads to the solution W\ = ctan(cry) and from the "bosonic" Riccati equation 



Wl + W[ + c\ + c 2 = , 



(17) 
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one will find cf (77) = 2c 2 sec 2 (c7y). Consequently, the Schrodinger-like equation 

-y" + cl(v)y = c 2 y (18) 

has solutions of the type y <x csec(c^), which is not of the localized soliton type, 
and therefore the approach leads to unphysical results. 

4 - Conclusion 

A cosmological KdV equation for adiabatic index solitons in the realm of barotropic 
FRW models has been introduced in this work. A simple recursive nonrelativistic 
SUSY scheme has been also applied to show that in the barotropic soliton profile 
one can embed a discrete sequence of adiabatic indices. Each eigenmode of the 
soliton corresponds to one of the embedded adiabatic indices. 
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